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Abstract. We present a useful method for the construction of cosmologi- 
cal models by solving the differential equations arising from calculating the 
kinematical invariants (shear, rotation, expansion and acceleration) of an ob- 
server field in proper time description. As an application of our method we 
present two generalizations of the Godel spacetime that follow naturally from 
g ; our approach. 

1. Introduction 

The construction of viable cosmological and astrophysical models often requires 
particular restrictions on the kinematic properties of those models. For example, 
O" 1 parallax-free models must necessarily be shear- free [8]. 

Here the cosmological observer field is analyzed in proper time description to 
obtain expressions for the metric that depend explicitly on the kinematical in- 
variants. This approach leads to a useful toolkit for constructing cosmological 
models with given kinematical properties. In addition, the analysis given here 
leads to a deeper understanding of the kinematical invariants and the relations 
between them. Of course, restrictions regarding the kinematical quantities also 
give rise to a limited range of possible matter models - it is often not possible 
to give energy-momentum tensors representing simple matter like a perfect fluid 



O 



X 



with a particular combination of rotation, shear, expansion and acceleration. Re- 
call for example the well-known shear-free fluid conjecture which has been proven 
for a number of special cases (cf. [7j). 

As an application of our method, we will construct a spacetime that general- 
izes the well-known Godel metric [1] which represents a dust model with negative 
cosmological constant that has non- vanishing rotation, but vanishing shear, accel- 
eration and expansion. It can be shown that this model contains closed timelike 
curves. Following [3], p. 170 those properties suggest that the Godel spacetime 
is not very physical. In 1952, Godel [2] made some remarks about the necessity 
of constructing a solution with non-vanishing expansion, but without giving an 



explicit metric. In this paper, we will construct one generalization of the Godel 
spacetime for which the acceleration vanishes as well as another model that is 
parallax-free. 

The authors wish to thank Professor Yuri Obukhov for many valuable suggestions 
while preparing the paper. 

2. Preliminaries 

We consider an (TV + l)-dimensional smooth Lorentzian manifold (M, g) of sig- 
nature (+, — ). Given coordinates ( X >j • • • • X ), by di we denote the partial 
derivative with respect to the Z th coordinate and Vj stands for the covariant de- 
rivative in direction of the j th coordinate vector. Latin indices will take the values 
0, . . . , N, whereas greek indices will range from 1 to N. In the following we fix a 
timelike unit vector field X 1 , the cosmological observer field. The usual decom- 
position of the covariant derivative of Xi into irreducible parts (see for example 
[I]) reads 

(1) VfcXj = u ik + a ik + —QP ik — XiX k 

with the antisymmetric part u)i k (rotation), the symmetric traceless part oi k 
(shear) and the trace G (expansion) itself. In detail, the parts read 

(2) uj ik = V[iXk] — X[iXk\, 

(3) o ik = V(iJV fc ) — X(iX k ) — — QP{ k , 
and 

(4) G = V a X a . 

Here the brackets (parentheses) denote the antisymmetric (symmetric) parts and 

(5) Pi k = gik — XiX k 

is the projection tensor on the TV-dimensional subspace perpendicular to X % . The 
acceleration is given by 

(6) X t = X a V a X t . 
In addition we have 

(7) X a X a = 0, a ai X a = u ai X a = P m X a = 0. 
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3. Proper time description of the observer field 



We adopt comoving coordinates with respect to the observer field such that X 1 = 
8q; this is also known as the proper time description. From X 1 being a timelike 
unit vector field we infer 1 = X a X a = g ab X a X b = goo- The first coordinate 
denotes the proper time of the cosmological observer; a fact which we will from 
now on emphasize by denoting x° with t. 

In proper time description the components of the metric can be expressed by 
components of the observer field - this is crucial to our analysis and we have: 



(8) 



Xi — g a iX a — g 0i 



and 



(9) 



Pik — Qik — goigok- 



For the covariant derivative one has 



(10) 



VkXi — dkXi — T a ik X a 
= d k Xi — T a ik g 0a 

= 9 k g 0i - -g ab g 0a (d k gib + dig kb - d b g ik ) 

= d k g 0i - ^5 b {d k g ib + dig kb - d b g ik ) 
1 

= o ( d k9Qi - dig ok + d g ik ) . 



For the acceleration, one has in particular: 



(11) 



Xi — X a V a Xi 



= V X, 

— dog 0i — -dig 00 

= d g 0i . 
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In the same manner, we obtain similar expressions for the other kinematical 
quantities in proper time description: 

(12) 6 = l -g %a d o9ia 

= ^d (logdet (g ik )) , 

(13) co ik = -(g kdo9oi- 9oid 9ok- d k goi + dig k), 
a ik — 2^°9 ik ~ 9okdo9oi — goidogok) — j^QPik 

(14) = \d P ik -^QP lk . 

It can also be easily seen from (J7|) that the kinematical quantities are characterized 
by their spatial components: 

(15) X = 0, uj 0i = a 0i = P 0i = 0. 

4. Models with given kinematical invariants 

Equation f Tl4l can be solved for the projection tensor P ik = g ik — goigo k , which 
leads to 

(16) P ik (t, x<) = x r ) + P lk (0, x 1 )) S 2 (t, x 1 ). 
with 

t 

1 



(17) S(t,x^)=exp | - / Q(T,xf)dr 



and 



t 



o 



Denoting derivation with respect to t by a dot, we have 

(19) e = jv| 

with S(0, x 1 ) = 1; 5(t,a; 7 ) is the so-called scale parameter. 



We interpret this equation as an evolution equation for the spatial components 
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of the metric for a given "start metric" ^ a/ g(0,x 7 ) and a given observer field 
9oa{t, xl )i expansion 6(£,x 7 ) and shear a aj g(£,x 7 ) that drive the evolution: 

(20) g al3 (t,x' Y ) = g 0a (t,x' Y )g 0l3 (t,x' Y ) + 

S 2 (t, x 7 ) (g af3 (0, x 7 ) - 90a (0, x 7 )^ 0/3 (0, x 7 )) + 
S 2 m aP (t,x<). 

Remark: The initial value u t = 0" is not a specific coordinate value like e.g. 
a singularity. In fact, the construction is not valid for singular start metrics. 
One should regard <7 Q/ g(0,x 7 ) as the spatial metric at some arbitrary time in the 
history of the observer. 

5. Some special cases 

Whereas expansion, shear and the observer field itself enter the evolution equa- 
tion directly, the rotation and acceleration impose additional constraints on the 
observer field. 

Models with vanishing acceleration. In the case of a vanishing acceleration, 
we have d g 0a = 0, which implies go a (t, x 7 ) = #oa(0, a; 7 ), and f[2"U]) reduces to 

(21) g a p{t,x*) = S 2 (t,x^)g af3 (0,x^ + 

(l-S 2 (t,tf))g 0a (0,x' r )gop(0,x' r ) + 
,S 2 (t)£ a/3 (t,x 7 ). 

Irrotational models. The rotation is calculated from the go a components of 
the metric alone. The equation for an observer field with vanishing rotation 
can be solved by setting g 0a = — /i(£)<9q,0(x 7 ). If the model also has vanishing 
acceleration, h(t) is constant and the observer field is a gradient. 

Parallax-free models. Spacetimes with an observer field parallel to a conformal 
vector field are important since such models are precisely the parallax-free models. 
As was shown in [8], this condition holds if and only if X 1 is shear-free and the 
exterior derivative of Xi — j^Xi vanishes. The last condition can be met by 
assuming 

(22) X t = ^X> - dj 
for some function /. 
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Integrating ( |22|) . and substituting into ( !20l) with a a/ 3 = 0, we have 

(23) g a/ 3(t,x r ) = S 2 (t,x'<)g a/3 (0,x'<) + 

S 2 (t,x~<)F a (t,x~ 1 )Fp(t,x~<) - 

S\t, x<) (F a (t, x~<)g 0/3 (0, x 7 ) + F p (t, x^)g 0a (0, x 7 )) 



with the functions 
(24) 



F n \t. X~* 



d a f(r, x 7 ) 
S(t, xt) 



dr. 



Toolkit. Figured] can be used as a toolkit to construct spacetimes with specific 
kinematic properties. We have only listed shear-free models since a non-vanishing 
shear can be easily included via the functions Also, we have notationally 

suppressed the dependence on the spatial coordinates x 1 . 



6. Generalizations of the Godel spacetime 

We like to generalize the Godel spacetime (R 4 , g) equipped with the usual observer 
field to a shear-free model with non-vanishing expansion. Written in canonical 
(cartesian) coordinates (t,x,y,z), the Godel metric reads: 



(25) 



(9 



ik 



( 1 






-1 



0\ 




U!()X 



V o 












V 



In this coordinates, the observer field is already of the desired form, X 1 = S l . We 
take as a start metric 



f-l 



(26) 



(^(O,^ 7 )) 



V 



o 








l p 2y/2uiQX 



V 



and as the observer field we choose 

(27) (g a(t,x~<)) = (o f(t)e^ o). 

The time-dependent function f(t) accounts for a possible no n- vanishing acceler- 
ation. 
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Spatial metric g a p(t) = gapit^x 1 ) 



General 



S 2 (t)g Q p(0) + g 0a (t)g p(t) - S 2 (t)g 0a (0)g p(0) 



S 2 (t)g a p(0) + (1-S 2 (t))g 0a (0)g p(0) 



6 = 



g a p(o) + go<x(t)goa(i) - go a (o)g 0a (o) 



= 0,X; = o 



gap{o) 



uj ik = 



S 2 (t)g a0 (O) + (h 2 (t)-h 2 (O)S 2 (t))d, 



Parallax-free 



S 2 (t)g af s(0) + S 2 (t) (F a (t)Fp(t) - F a (t)g op (0) - F p (t) g0a (0)) 



Figure 1. Table of shear-free cosmological models 



In order to evolve this metric in a shear-free manner we use equation fl2"U]) with 
°"a/3 = 0. To simplify things, we assume that the scale parameter S depends only 
on the cosmological time parameter t. In this way we arrive at 
(28) 

(-S 2 (t) \ 



(#a/3(M 7 )) 



MS 2 (t)(l-2/ 2 (0)) + 2/ 2 (t)]e 







V 











-S 2 (t)J 



We like to model the matter content of this spacetime by a fluid (not necessarily 
perfect). Thus, the energy-momentum tensor takes the form 



(29) 



T ik = -pgik + (e + P)XiX k + 2q (i X k) + 7v ik , 
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with energy density e, isotropic pressure p, heat flow and anisotropic pressure 

The heat flow and anisotropic pressure have to satisfy the conditions q a X a = 0, 
7r a iX a = and n a a = 0. The first condition directly implies qo = 0, whereas the 
second condition fixes g 7 by inspection of the Einstein equation Rik — \Rg%k — 
^■9ik — KTik after solving for e and p: 



(30) * = M2/t-i) i|( 3/( ' )s(i) ' s( ' >/( ' 1 )' 

(31) g 2 = ~ \g 2(t) {S(t)S(t)-S(t)), 

(32) g 3 = 0. 

Models with vanishing acceleration. First we deal with the case of vanishing 
acceleration and set f(t) = 1. The evolved metric then reads 

/ 1 e^ 0X \ 

-S 2 {t) 



(33) (g, 



e V2u x o l e 2sfiu x ( 2 -S 2 (t)) 
\ S 2 {t)J 

and the condition 7r a a = enforces 

(34) 2u 2 (1 - 5 2 (t)) - S(t)S(t) + 2(5(t)) 2 = 0. 
This equation can be reduced to the separable first order equation 

(35) (S(t)) 2 = ^o{^~^S 4 (t) + 2S 2 (t)-l) 
which implies 

(36) S(t) = 2u 2 S{t) ((a - 1) S 2 {t) + l) . 

Substituting ( J35l) and ( J36l) into the differential equation (!34l) one sees that it is 
identically satisfied for any value of the constant a. Furthermore, the equation 
can be solved in terms of elliptic integrals (the solution is one of the Jacobi elliptic 
functions). 

However, evaluating fl33|) at t — one sees that 9au; 2 = 9S' 2 (0) = O 2 (0) and we 
must have a > to obtain real valued solutions. For a = = O(0) one has 
S(t) = 1 and we recover the Godel metric. The phase portrait of the differential 
equation for different values of 9(0) is shown in figure EJ If < |6(0)| < 3|o>o|, the 
solutions are periodic, whereas for 3\u \ < |6(0)| the spacetime expands/collapses 



indefinitely. The metric does not become singular at any point since S(t) does 
not vanish for any initial value. 

Furthermore, the model exhibits an interesting causal structure that we will dis- 
cuss in the following section. 




Parallax-free models. We will now choose f(t) such that the model is parallax- 
free. Calculating the exterior derivative of Xi — %Xi, we find that the observer 
field of our already shear-free model (12 81) is parallel to a conformal vector field if 
and only if the following set of differential equations is satisfied: 

(37) f(t)S(t)-S(t)f(t) = 0, 

(38) f(t)S(t)-S(t)f(t) = 0. 

Together with the condition of a traceless anisotropic pressure we finally have 

S(t) = e5 0t = f{t) with the constant expansion O. 
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The resulting metric takes the form 



(39) 



(9ik) = 



( 



\ 



1 






_ e §et 








1 p 2V2ui x+%Bt 
2 e 3 









\ 







-e>t 



which is precisely the Godel metric for = 0. 

We compute some scalar invariants of this solution: Rotation scalar and the norm 
of the acceleration are given by 



ab 



et 
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-e 2 . 



(40) u ab u" 

(41) X a X a 

The heat flow and anisotropic pressure are non-zero: 
(42) 



q a q = K ah -K = - 



8lu 2 Q 2 



:0t 



9 K 2 



Finally, the energy density and pressure are given by: 

(43) e = i- (3 Wo 2 e-S et - 2 - 3A) , 

(44) p 



1 

3k 



3uj 2 e-i et + e 2 + 3A 



A particular interesting case is A = —^y, where we have a stiff matter model 
(e = p — ^e - ^ *) with non-vanishing anisotropic pressure and heat flow that 
is asymptotically a lambda vacuum for large cosmological times. In general, 



one must have A < 



e 2 



to ensure that the cosmological observer measures an 



everywhere non-negative energy density T ab X a X b = e. 



7. Closed causal curves 

Any function on a smooth closed curve x l (s) is periodic with respect to the curve 
parameter s. If the curve is contained in a chart of a Lorentzian manifold, this 
holds in particular for the coordinate function x°(s) and there exist maximal and 
minimal points of x°(s) for which ^ = holds. Let x l (s) be such a closed curve 
that is also causal, i.e. 

(45) *-5"sr £ 
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holds for all s. If So is chosen such that ^ = holds (which is always possible 
by the arguments above) then 

dx a dx 13 



, , , dx 1 dx k 

(46) g ik 



ds ds 



^ al3 ds ds 

s=so 



S = SQ 



is negative if the matrix consisting of the g a p is negative definite - which therefore 
is a sufficient condition for the absence of closed causal curves. 
Following basic linear algebra a real square matrix A = [a a p] is negative definite 
if — A is positive definite. According to the Sylvester criterion, this holds if and 
only if the determinants of all leading principal minors — A^ (beginning with the 
element —an) are positive. 

Examining the generalized Godel metric with vanishing acceleration ( 1331) . one 
immediately sees that this condition is met for S 2 (t) > 2. Thus, if the initial 
expansion compared to the initial rotation is large enough such that the solution of 
( |34|) takes large enough values, the cosmological observer in the acceleration-free 
expanding generalized Godel spacetime "escapes" its causality violating epoch. If 
H | Wo | < | © (0) | < 3|o;o|, the observer periodically enters the chronology violating 
region of spacetime. If 3 \uq\ < |@(0)|, this region is eventually avoided. Note 
that our argument does not show that the causality violating region is as large 
as {S 2 (i) < 2}. In particular, y| |o> | < | O (0) | may not be the best estimate for 
the initial values that give rise to causally well-behaved regions of this particular 
spacetime. 



8. Conclusion 

The construction of cosmological models with given kinematical invariants is an 
important task. On the one hand, it is an important problem in mathematical 
cosmology to examine possible connections between the kinematics of a cosmo- 
logical model and its geometric properties like singularities and causality. On the 
other hand, astrophysical observations put constraints to the possible kinematical 
quantities of our universe. With the description of spacetime metrics presented 
here, these constraints can be directly accounted for. 

The usefulness of our method is clearly demonstrated by the given examples. 

The parallax-free generalization of the Godel metric presented here seems to be 

more natural in terms of symmetry since we have just relaxed the property of a 

Killing observer field to an observer field parallel to a conformal vector field. The 
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model with vanishing acceleration on the other hand has an interesting geometric 
structure with causality violating regions and periodic solutions. 
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